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Abstract 



. . We report the first observation of the decay B^ — > XcqK'^ using 21.3 fb^^ 

Q I of data collected by the Belle detector at the T(4S) resonance. The prelimi- 

^^ I nary result for the branching fraction is B{B^ -^ XcqK^) = (8.0_2'4 ± 1.0 it 

r~^ ' 1.1) X 10^^ where the first error is statistical, the second is systematic, and 

^D ■ the third conies from the uncertainty in the XcO ~^ tt+tt^ branching fraction. 
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I. INTRODUCTION 

The production rate of quarkonium states in various high energy physics processes can 
provide valuable insight not only into the interactions between a heavy quark and antiquark, 
but also into the elementary processes with the production of a QQ pair. The spin-1 5"- 
wave resonances, such as the J /ip of the charmonium family, are of special experimental 
significance because of their importance for studies of CP violation. The P states are also 
important in their own right because they probe a qualitatively different aspect of the QQ 
production process. While the S states probe only the production of color-singlet QQ pairs 
at small separation distances, the P states can also probe the production of color-octet QQ 
pairs. Color-singlet production of Xco iii B decays vanishes in the factorization approximation 
as a consequence of spin-parity conservation. However, the color-octet mechanism allows for 
the production of the Xco -P-wave 0"*"+ state via the emission of a soft gluon [0,|| . 

The only recent search for the B -^ XcoK decay was reported by CLEO ||^. Using 
9.66 X 10^ BB pairs they set 90% C.L. upper limits of: B{B+ -^ XcoK+) < 4.8 x 10"^ and 
B{B^ -^ XcqK^) < 5.0 X 10"'^. This analysis uses a data sample collected with the Belle 
detector [Q at KEKB asymmetric energy e~^e~ collider [^. It consists of 21.3 fb~^ taken at 
the T(45') (which corresponds to about 22.8 million produced BB pairs) and 2.3 fb^^ taken 
60 MeV below for continuum studies. 

The detailed analysis of the final state presented in this paper is part of the general 
analysis of i?-meson three-body decays reported in ref. [Q. 

The inclusion of charge conjugate states is implicit throughout this report unless explic- 
itly stated otherwise. 

II. EVENT SELECTION 

Charged tracks are required to satisfy a set of track quality cuts based on the average hit 
residual and impact parameters in both r-0 and r-z planes. We require that the transverse 
momentum of the track be greater than 100 MeV/c to reduce low momentum combinatorial 
background. For more details, see ref. [Q. 

Hadron identification is accomplished using the responses of the ACC and the TOF and 
dE/dx measurements in the CDC. The information from these three subsystems is combined 
in a single number using the likelihood method: 

c{h) = c^^'^ih) X /:™^(/i) X c^'^'^ih), 

where h stands for the hadron type (tt, K, p). Charged particles are identified as K^s or vr's 
by cutting on the likelihood ratio (PID): 



The likelihood ratio for kaon candidates is required to be greater than 0.5 and for pion 
candidates less than 0.9. 

All charged tracks are also required to satisfy an electron veto requirement that demands 
that the electron likelihood is less than 0.95. In addition, all charged kaon candidates are 
required to satisfy a proton veto: 

For more details see rf. [§ and references therein. 

We reconstruct Xco candidates in the Xco — ^ vr+vr" and Xco — *■ K'^K' decay modes. The 
B^ -^ XcqK^ candidate events are identified by means of the beam-constrained mass Mbc 
and the energy difference AE: 



Mbc = ^s/A - P^^ ; AE = E*^ - v^/2, 



where E^ and P^ are the energy and three-momentum of the B candidate in the T(4S') rest 
frame and ^/s is the total energy. Subsequently, we refer to the "5 signal region," which is 
defined as: 

5.272 < Mbc < 5.289 GeV/c^; \AE\ < 40 MeV. 

III. BACKGROUND SUPPRESSION 

To suppress the combinatorial background which is dominated by the two-jet-like 
e~^e~ —>■ qq continuum process, we use variables that characterize the event topology. We 
require | cos(6't/it-)| < 0.80 where Oxhr is the angle between the thrust axis of the B candidate 
and that of the rest of the event. This eliminates 83% of the continuum background and 
retains 79% of the signal events. Following the analysis of ref. Q, we define a Fisher discrim- 
inant JF that includes the angle of the thrust axis of the entire event, the production angle of 
the B meson candidate and nine "Virtual Calorimeter" parameters that characterize the 
momentum flow in the event relative to the B candidate thrust axis. When the requirement 
on the Fisher discriminant variable JF > 0.5 is imposed, 79% of the continuum background 
is rejected with about 74% efficiency for the signal. In the case of three charged kaon final 
states, the continuum background is much smaller and a looser requirement JF > can be 
imposed. This rejects 53% of the continuum background with about 89% efficiency for the 
signal. 

In the A'+TT+TT^ final state the background from B generic decays is dominated by the 
decays of type B^ -^ [A'+7r^]7r+ where [A'+tt^] denotes a resonance state which can decay 
into A'+TT" such as D^, A'*°(892) etc. To suppress this type of background, we require 
the invariant mass of the A'+vr" system be greater than 2.0 GeV/c^. For the K^K^K~ 
final state we require the invariant mass for both K^K~ combinations be greater than 2.0 
GeV/c^ to suppress the possible background from rare B decays such as B^ -^ (pK^. 



IV. RESULTS OF THE ANALYSIS 

A. Signal yield extraction 

First, we select all K+tt+tt^ (K^K^K^) combinations from the B signal region that sat- 
isfy the selection criteria described above and have the vt+tt^ {K'^K~) invariant mass in the 
range 3.2 < M{h~^h~) < 3.8 GeV/c^. The resulting vr+vr^ and K^K^ invariant mass spec- 
tra are shown in Figs. |I|a and |T]b, respectively. Since in the case of the K^K^K~ final state 
there are two same charge kaons, we distinguish the K^K~ combination with the smaller, 
M[K^K~)min-, and larger, M{K^K~)maxi invariant masses. Only the larger combination 
is plotted. Two prominent peaks can be seen in Fig. |l]a. The peak around 3.69 GeV/c^ 
corresponds to the ip{2S) meson from the decay mode B^ — >■ ip{2S)K^, ip{2S) -^ /i"*"/!" 
with muons misidentified as pions. It is slightly shifted to higher mass values due to the 
wrong mass assignment for ip{2S) daughter particles. The peak just above 3.4 GeV/c^ is 
identified as the Xco meson. 

The TT+vr" and K^K~ spectra are fitted to the sum of a Breit-Wigner function convolved 
with a Gaussian for the signal and a zeroth-order polynomial for the background. The width 
of the Gaussian is fixed at 10.8 MeV/c^ as determined from a fit to the J /ip peak in the 
/i"^/i~ invariant mass spectrum. The width of the Breit-Wigner function is fixed at the 
world average value for the Xco IK- The results of the fit to the tt+tt^ and K^K~ invariant 
mass spectra are given in Table |I[ The peak position in the K^K^ spectrum. Fig. |p3, is 
found to be somewhat shifted to lower mass values. Although this shift is consistent with 
a statistical fluctuation, we note that according to the results of our general analysis of the 
three body B^ -^ K^K^K~ decay, we find a significant signal in the mass sidebands of 
the Xco ~^ K^K~ signal 0. As a result, the K^K~ invariant mass distribution in the Xco 
region could be distorted by the effects of interference with an amplitude not related to the 
B^ -^ XcqK^. a Monte Carlo simulation study indicates that this effect could shift the Xco 
peak position by as much as 15 MeV/c^ and could have as much as a 100% effect on the 
observed amphtude of the Xco signal. Because of this uncertainty we base our branching 
fraction measurement for the B^ —>■ XcqK^ decay on the Xco ~^ tt+tt^ decay mode only. 

TABLE I. Results of the fit to the vr+vr^ and K^K^ invariant mass spectra. 



XcO submode 


Efficiency, % 


Peak, GeV/c2 


Yield, events 


Significance, a 


XcO -^ vr+vr- 
Xco -^ K+K- 


21.5 
13.7 


3.408 ± 0.006 
3.390 ± 0.009 


15.511^ 
7.7tl:? 


4.8 
3.2 



The XcO candidates are then selected by requiring \M{h^h^) — 3.415| < 0.050 GeV/c^, 
where h^h' denotes both vr+vr" or K^K~ combinations. This corresponds to about a 2.6 a 
cut. The resulting two-dimensional /S.E versus Mbc plots as well as the projections on the 
/S.E and Mbc axes are shown in Figs. |^ and ^ for the K^n^n' and K^K^ K~ final states, 
respectively. Clear signals are apparent in both figures. 



B. Branching fraction calculation 

To determine branching fractions, we normalize our results to the observed B^ —>■ D^7t~^, 
D^ —>■ K^TT~ signal. Although this introduces a 9.7% systematic error because of the 
uncertainty in the B^ -^ D^n^ branching fraction, it removes systematic effects in the 
particle identification efficiency, charged track reconstruction efficiency and the systematic 
uncertainty due to the cuts on event shape variables. We calculate the branching fraction 
for the B meson decay to a final state / via the relation: 

B{B ^ f) = B{B+ -^ D\+) X i3(D° ^ K+n-)^ x ^, (1) 

where Nf and Non are the numbers of observed events for the particular final state / and for 
the reference process respectively, Ef and eon are corresponding reconstruction efficiencies 
determined from the Monte Carlo simulation. 

As a cross-check, we use the decay mode B^ -^ J /ipK^ followed by the J /ip -^ fi^^~ 
decay. In this case there are two muons in the final state instead of two pions. To avoid 
additional systematic uncertainty in the muon identification efficiency, we do not use muon 
identification information for J/ip reconstruction. Instead, we apply the same pion-kaon sep- 
aration cut for muons from J/ip as for pions from Xco- The feed-down from the J/ip — > e^e~ 
submode is found to be negligible (less than 0.5%) after the application of the electron veto 
requirement. 

The signal yields for the B^ —>■ D^n^ and B^ —>■ J/tpK^ processes are extracted from 
fits to the AE distributions shown in Fig. §. The yields are found to be Nd.„ = 858 ±31 and 
Nj/^ = 323 ± 19 respectively. Using the corresponding reconstruction efficiencies of 23.5% 
and 29.4% we find the ratio, R, of branching fractions to be 

where only the statistical error is quoted. This result is in good agreement with the value 
calculated from PDG |§ data: 0.189 ± 0.026. 

To calculate the branching fraction for the B~^ — > XcqK^ decay mode, we use the signal 
yield determined from the fit to the h~^h~ invariant mass spectra to take into account the 
possible contribution from the non-resonant i?+ —>■ K~^h^h~ decays that would produce 
signal-like distributions in both AE and Mbc- Combining all the relevant numbers from 
Table | and the intermediate branching fractions from PDG [Q, we find the branching 
fraction for the B^ — > XcqK^ decay mode to be: 

B{B+ -^ XcoK+) = {S.Otli ± 1.0 ± 1.1) X 10-^ ixco ^ vr+vr- mode); 

B{B+ -^ XcqK+) = (5.3i|^ ± 0.7 ± 0.8) x 10"^ {xco ^ K+R- mode), 

where the first error is statistical, the second is systematic, and the third is due to the 
uncertainty in the XcO ~* h~^h~ branching fractions. The systematic error consists of the 
uncertainty in the B^ -^ D^n^ and D^ -^ K^tt^ branching fractions (9.7%), the uncertainty 
in the background parameterization in the fit to the h^h' spectra (7.8% for tt^tc" and 8.5% 
for K^K~) and the uncertainty in the background and signal parameterization in the fit to 



the AE distribution for B^ — > D^tt"*" signal (2.3%). As mentioned above, our qualitative 
estimates indicate that the model dependent error in the K^K^K~ final state could be very 
large. We do not present the quantitative estimation here because of the large uncertainty. 
We also calculate the ratio of the branching fractions for the B^ —>■ XcoK^ and 
B^ -^ J l^ipK^ decays using only the Xco ^ tt+tt^ submode: 

where the first error is statistical and the second comes from the uncertainties in the 
J /ip -^ fi'^fi~ and Xco -^ vr+vr" branching fractions. 

C. Cross-check with XcO —>■ K*^K^tt^ mode 

As an additional cross-check, we reconstruct the Xco nieson in the multi-body 
K^K^Ti^Ti^ final state. Although, in general, the branching fractions for Xco decays to 
multi-body final states are significantly higher than those for tt+tt" and K^K~ , these modes 
suffer from a much larger combinatorial background. This background is somewhat reduced 
in the K~^K~7r~^7r~ final state because of the presence of two charged kaons. We further 
reduce the combinatorial background by requiring that at least one Ktt pair be consistent 
with K*{892) -^ Kn: \M{Ktt) - 0.896| < 0.050 GeV/c^. We also use a tighter cut on the 
Fisher discriminant variable JF > 0.8 to suppress the larger continuum background. 

The resulting two-dimensional AE versus Mbc plot as well as the AE and Mbc projec- 
tions are presented in Fig. ^ From the fit to the AE distribution 9.2i3;4 signal events are 
found. Using the reconstruction efficiency of 7.2% from Monte Carlo, we obtain a branching 
fraction that is in agreement with that determined for the vr+vr" mode, namely 

l3iB+ ^ XcoK+) = (7.011° ± 2.3) x 10"^ (xco -^ K*^K~n+ mode), 

where the first error is statistical and the second comes from the uncertainty in the 
XcO -~^ K*^K~7T^ branching fraction. 



V. CONCLUSION 

We report the first observation of the B^ -^ XcqK^ decay mode. The statistical 
significance of the signal is more than 4o". The preliminary branching fraction result is 
B{B+ -^ XcoK-^) = (S.OllI ± 1.0 ± 1.1) X 10-^ which is comparable to that for the 
5+ —>■ J/tpK'^ decay. This provides evidence for a significant nonfactorizable contribution 
in B to charmonium decay processes. 

The branching fraction measured by Belle is slightly higher than the CLEO published 
upper limit. However, given the uncertainties on the result, it would be premature to claim 
that they disagree. 



ACKNOWLEDGEMENT 

We wish to thank the KEKB accelerator group for the excellent operation of the KEKB 
accelerator. We acknowledge support from the Ministry of Education, Culture, Sports, 
Science, and Technology of Japan and the Japan Society for the Promotion of Science; 
the Australian Research Council and the Australian Department of Industry, Science and 
Resources; the Department of Science and Technology of India; the BK21 program of the 
Ministry of Education of Korea and the CHEP SRC program of the Korea Science and 
Engineering Foundation; the Polish State Committee for Scientific Research under contract 
No.2P03B 17017; the Ministry of Science and Technology of Russian Federation; the National 
Science Council and the Ministry of Education of Taiwan; the Japan- Taiwan Cooperative 
Program of the Interchange Association; and the U.S. Department of Energy. 



REFERENCES 

[1] G. Bodwin et al, Phys. Rev. D 46, 3703 (1992). 

[2] M. Beneke et a/., Phys. Rev. D 59, 054003 (1999). 

[3] K.W. Edwards et al. (CLEO Collaboration), Phys. Rev. Lett. 86, 30 (2001). 

[4] K. Abe et al. (Belle Collaboration), KEK Progress Report 2000-4 (2000), to be pubhshed 

in Nucl. Inst, and Meth. A. 
[5] KEKB B Factory Design Report, KEK Report 95-7 (1995), unpublished; Y. Funakoshi 

et al, Proc. 2000 European Particle Accelerator Conference, Vienna (2000). 
[6] K. Abe et al. (Belle Collaboration), BELLE-CONF-0114; Submitted as a contribution 

paper to LP2001. 
[7] D.M. Asner et al. (CLEO Collaboration), Phys. Rev. D 53, 1039 (1996). 
[8] D.E. Groom et al. (Particle Data Group), Eur. Phys. J. C 15, 1 (2000). 



10 



12 



10 



8 



> 
in 



o 6 - 



(O 



o 
o 




3.2 3.3 3.4 3.5 3.6 3.7 3.8 

M(7i;V) (GeV/c^) 



*y 3 



o 

lO 



w 2 



o 
o 




3.2 3.3 3.4 3.5 3.6 

M(K^K) (GeV/c^) 

FIG. 1. The vr^TT^ (a) and the K^K^ (b) invariant mass spectra for B candidates from the 
signal region. The histograms are data and the curves are fits to the data. 
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FIG. 2. The AE and M^c distributions for S+ 
the two-dimensional plot shows the signal region. 
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FIG. 3. The /S.E and Af(,c distributions for B^ — > XcqK'^ , XcO -^ K^K candidates. The box 
in the two-dimensional plot shows the signal region. 



13 



.750 




0.1 0.2 

AE, GeV 



60 






§ 45 



30 



15 



+^_ ++Uj 




%.2 -0.1 0.1 

AE, GeV 

FIG. 4. The AE distributions for reference processes: a) - B^ -^ D°7r+, D" 
h) - B+ ^ J/^K+, J/^P -^ n+fi-. 



K+TT^; 



14 



0.2 



0.1 



> 

a 

<l 



-0.1 



-0.2 



.■■ .~- ■ . 

- "■ ■ 

-■ " ■ ■ " 

■ ■ 

■ ■ 

■ ■ ■ ^ ■ ■ ■ 

- ' ■■■ .* ■ ■.■'■ 

■ . . .■ ■ 

"■ ■ ■ ■ 

■ ■ 

■ *.-."■ 

1 ■ ■ 
■ ■ _ ■ ■ 

■■ ■■ :■. ■■ ■ 
■■' -^ . ■■ .■— .■ 

■ i ■ ". ■■ 

.: ■': ■ .■■> . ■■■ . y 
■■. ■■■■.; ■ ' .'■■.■.: 


■ " 
■ 








■ ■ 

■ ■ 
...■■ ■ . 

■i-, ■-, 





5.2 



5.225 



5.25 
Mbc (GeV/c^) 



5.275 



5.3 





Mbc (GeV/c^) 

FIG. 5. The l\E and M^c distributions for 5+ ^ XcO-?^+, XcO 
box in the two dimensional pfot shows the signal region. 
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